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We investigate optimal separable approximations (decom- 
positions) of states £1 of bipartite quantum systems A and 
B of arbitrary dimensions M x A'^ following the lines of Ref. 
[M. Lewenstein and A. Sanpera, Phys. Rev. Lett. 80, 2261 
(1998)]. Such approximations allow to represent in an opti- 
mal way any density operator as a sum of a separable state 
and an entangled state of a certain form. For two qubit sys- 
tems (Af — N = 2) the best separable approximation has a 
form of a mixture of a separable state and a projector onto a 
pure entangled state. We formulate necessary condition that 
the pure state in the best separable approximation is not to 
be maximally entangled. We demonstrate that the weight of 
the entangled state in the best separable approximation in 
arbitrary dimensions provides a good entanglement measure. 
We prove in general for arbitrary M and N that the best sep- 
arable approximation corresponds to a mixture of separable 
and entangled state both of each are unique. We develop also 
a theory of optimal separable approximations for states with 
positive partial transpose (PPT states). Such approximations 
allow to decompose any density operator with positive partial 
transpose as a sum of separable state and an entangled PPT 
state. We discuss procedures of constructing such decompo- 
sitions. 

03.67.Hk, 03.65.Bz, 03.67.-a, 89.70.-hc 



I. INTRODUCTION 

The problem of characterization of entangled states of 
composite quantum systems is one of the fundamental 
open problems of quantum theory. Entanglement is one 
of the quantum properties which make quantum mechan- 
ics so fascinating: it leads to famous apparent paradoxes 
[Q,^, and it is of great importance for applications in 
quantum communication and information processing . 

In the case of the pure states it is easy to check whether 
a given state is, or is not entangled. So far, the answer to 
this question when applied to quantum mixtures is not 
known in general. The definition (introduced by Werner 
Q ) says that a state (in general a mixed state) is entan- 
gled when it is not separable. Separable states defined 
on a Hilbert space Hab — 'Ha ®'Hb are those that can 
be as a convex combination of projections onto product 
states 
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In finite dimensional spaces, the number of terms in 
the sum can be restricted to < dim(HAs)^ (in another 



words, when the density matrix is separable, then it can 
be represented in the above form with K terms, where 
K is not larger than the dimension of the space of linear 
operators acting in TLab, see (H). 

Several necessary conditions for separability are 
known: Werner's condition based on the mean value 
of the, so called, flipping operator [Q, Horodeckis cri- 
terium based on a-entropy inequalities Q , and many oth- 
ers 0. Perhaps, the most important necessary criterium 
has been formulated by Peres who has demonstrated 
that the partial transpose g^^ of any separable matrix 
Q defined as {m, iJL\Q^^\n,v) = (n, /i|£)|m, i/) for any fixed 
orthonormal product basis | 71, v) = |cyj)^ ® 1^1^) B must be 
positively defined. In the following we will call states with 
positive partial transpose PPT states. Physical mean- 
ing of the PPT property is for PPT state time reversal 
operation in one subsystem (either Alice's or Bob's) is 
physically sound 

It is worth stressing that the problem of separability 
is directly related to the theory of positive maps on C*- 
algebras [|lO|,0 This has been established in Ref. 
in which it was shown in particular that for systems of 
low dimensions [M x N < 6) the PPT condition is also 
sufficient for separability. For systems of higher dimen- 
sions {M xN > 6) there exist entangled states having the 
PPT property. First examples of such were provided by 
means of the, so called, range separability criterion based 
on analysis of range of density matrix (see also |0 ) . 
Such states represent represent bound entanglement, i.e. 
cannot be distilled JTsf . 

In the recent Letter we have also looked at the range of 
the entangled density operators in order to formulate an 
algorithm of optimal decomposition of mixed states into 
the separable and inseparable part Our method 

of the best separable approximations (BSA) was based 
on subtracting projections on product vectors from a 
given density matrix in such a way that the remainder 
remained positively defined. This approach allowed to 
achieve a variety of vary strong results: optimal decom- 
positions with minimal number of terms in the form of 
mixtures and pseudo mixtures for 2x2 and 2x3 systems 
Q , separability criteria for 2 x systems [ and in gen- 
eral for M X N systems (with AI < N) |16| for density 
matrices of low ranks. In particular it was shown that: i) 
all PPT states of rank smaller than A^ are separable; ii) 
1. (1) for generic states such r{g) + r{g'^^ ) < MN -M -N + 2 

constructive separability criteria were derived that re- 
duce the problem to finding roots of some complex poly- 
nomials; iii) for 2 X A^ it was shown that for the states 
invariant under partial transpose with respect to the 2 
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dimensional subsystem, and those that are not "very dif- 
ferent" from their partial transpose are necessarily sep- 
arable. Very recently, these findings have allowed us to 
present general schemes of constructing non decompos- 
able entanglement witnesses (i.e. observables that have 
a positive mean value on all separable states, and have a 
negative mean value on a PPT entangled state ) and 
nondecomposible positive maps in arbitrary dimensions, 
that is the maps that cannot be decomposed into a sum of 
a completely positive map and another completely posi- 
tive map combined with the transposition [ p^ . It should 
be stressed that our approach goes beyond the methods of 
constructing examples of PPT entangled states and pos- 
itive maps based on the, so called, unextendible product 
bases [0,^ . More importantly, we were able to present 
methods of constructing optimal entanglement witnesses 
and optimal nondecomposible maps which provide very 
strong separability criteria |^ . In a series of importan- 
tant papers Englert and his collaborators have obtained a 
series of remarkable analytic results concerning the BSA 
decompositions for 2 x 2 systems |Q . These results give 
new deep insight into the fundamental problem of quan- 
tum correlations in 2 qubit systems. 

All of the above mentioned applications indicate that 
the method of BSA is very useful. The aim of this paper 
is to generalize and to complete results of the Refs. [ p^ . 
We present several results that characterize the BSA de- 
compositions in 2x2 and, in general in M x iV systems. 
Concerning the 2 qubit systems our results are comple- 
mentary to those of Ref. [|T]. The plan of the paper 
is as follows: In Section II we remind the reader some 
basic facts about the optimal and the best separable ap- 
proximations. In Section III (using also the results pre- 
sented in the Appendix) we demonstrate necessary condi- 
tion that for a two qubit systems {M = N = 2) the best 
separable approximation has a form of a mixture of a sep- 
arable state and a projector to an entangled state which 
is not maximally entangled. In Section IV we remind the 
reader the basic facts about entanglement measures; we 
prove here that the weight of the fully entangled state 
in the BSA decomposition of M x states provides a 
good entanglement measure. In Section V we prove that 
in general for arbitrary M and N the best separable ap- 
proximation corresponds to a mixture of separable and 
entangled state both of each are uniquely determined. Fi- 
nally, in Section VI we formulate the theory of optimal 
separable approximations for states with positive par- 
tial transpose (PPT states). Such approximations allow 
to represent any density operator with positive partial 
transpose as a sum of separable state and an entangled 
PPT state. Decompositions of this sort play essential 
role in the theory of nondecomposible positive maps [ p^ . 
We present and discuss efficient numerical procedures of 
construction of such decompositions. 



II. INTRODUCTION TO BSA 

Consider a state p acting on C*^ (E) . Such a state 
will be called a PPT state if its partial transpose sat- 
isfies p^^ > (or equivalently p^'^ > 0). Throughout 
this paper K{X), R{X), k{X), and r{X) denote the ker- 
nel, the range, the dimension of the kernel, and the rank 
of the operator X, respectively. By |e*) we will de- 
note the complex conjugated vector of |e) in the basis 
|0)a, • ■ • in which we perform the partial transposi- 
tion in the Ahce space; that is, if |e) = q;|0) + + . . . 
then |e*) = a*|0) + + . . .. Similar notation will be 
used for vectors in the Bob's space. 

In this section we give a short repetition of what we call 
optimal and the best separability approximations (OSA, 
and BSA respectively). Although the results below have 
been proven in Ref. jl^ , we repeat them here using the 
notation of the present work. The idea of BSA is that, 
because of the fact that set of separable states is compact, 
for any density matrix p there exist a "optimal" separa- 
ble matrix p* and "optimal" A > such that Ap* can 
be subtracted from p maintaining the positivity of the 
difference, p — Ap* > 0. This situation is characterized 
by the following theorem: 

Theorem 1 For any density matrix p (separable, or not) 
and for any (fixed) countable set V of product vectors 
belonging to the range of p, i.e. \ea, fa) G R{p), there 
exist h.{y) > and a separable matrix 

P:(\^) = ^A„P„, (2) 

a 

where Pa = \e,a, fa){e-a, fa\, while all A^ > 0, such that 
6p — p — Ap* > 0, and that p*g{V) provides the opti- 
mal separable approximation (OSA) to p since Tr(5p) 
is minimal or, equivalently A is maximal. There ex- 
ists also the best separable approximation p* for which 
A = maxyA(y). Obviously, A{V) < A{V') when V CV 

Remark 1 Quite generally one can define the best sepa- 
rable approximations ps of p by demanding that — ps|| 
is minimal with respect to some norm in the (Banach) 
space of operators. Here we minimize Tr(p — \ps) with 
respect to all ps such that p — Xps > 0. 

From this theorem it follows then that if any density 
matrix p is separable then A = 1. Caratheodory's the- 
orem implies then (see discussion in Ref. [1) that there 
exist a finite set of product vectors V C R{p) of cardinal- 
ity < r{p)^, for which the optimal separable approxima- 
tion to p, p* \y] is equal to the BSA and A = 1 also. The 
above theorems are also true for uncountable families of 
states V , and appropriate generalizations are discussed 
in Ref. 

In order to explain now how the procedure of construc- 
tion of the matrix p* actually works, we introduce two 
important concepts: 
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Definition 1 A non-negative parameter A is called 
maximal with respect to a (not necessarily normalized) 
density matrix p, and the projection operator P = IV'XV'I 
if p~AP > 0, and for every e>0, the matrix p—{ A +e)P 
is not positive definite. 

This means that A determines the maximal contribution 
of P that can be subtracted from p maintaining the non- 
negativity of the difference. Now we have the following 
important lemma: 

Lemma 1 A is maximal with respect to p and P = 
if: (a) if liP) ^ R{p) then A = 0, and (h) if 
IV') e R{p) then 



< A 



(3) 



Note that in the case (b) the expression on RHS of Eq. 
^ makes sense, since jV") G R{p)^ and therefore there 
exists 10) such that = p\4>), or equivalently that 
p~^|V') = !</>)■ Remarkerbly this Lemma has been used in 
a completely different context by E. Jaynes in his works 
on foundations of statistical mechanics |Q . 

Definition 2 A pair of non-negative (Ai,A2) is called 
maximal with respect to p and a pair of projection opera- 

tors Pi = |V'l)(^l|, P2 = |MV'2|, lfp-KiPi~K2P2 > 0, 

Ai is maximal with respect to p — A2P2 and to the pro- 
jector Pi, A2 is maximal with respect to p — AiPi and to 
the projector P2, and the sum Ai + A2 is maximal. 

The condition for the maximality of Ai + A2 is the given 
by the following lemma: 

Lemma 2 A pair (Ai, A2) is maximal with respect to p 
and a pair of projectors (Pi,P2) */•' 

• (a) if \ipi) , |V-'2) do not belong to R{p) then Ai — 
A2 = 0; 

• (h) if \%Iji) does not belong to R{p), while 1-02) G 
R{p) then Ai ^ 0, A^ = {^2\p-^\^2)-\- 

• (c) ifl^i), \ip2) e R{p) and (i/'i|p"^|'(/'2) = 0, then 
A, = (0,|p-i|i^,), 1 = 1,2; 

• (d) if IV-i), IV'2) e R{p) and 
(Vi|p-Vi),(V'2|p-^|V'2) > \{^Pi\p-'\^2)\ ¥^ then 

Ai = ((V'2|p-V2) - \{i^l\p~'\^/j2)\)/D, (4) 
A2 = ((V-ilp-Vi) - \{Hp-'\i^i)\)/D, (5) 



where 

D = (V'i|p^^|0i)(^2|p-i|02) - 1(011^-11^2)^; 

(e) finally, i/ |0i),|V'2) e R{p) and (0i|p"^|0i) > 
|(0i|p-i|V'2)| > {^2\p-^\iJ2), then Ai = 
{^Pl\p-^\^Pl)-\ A2 ^ 0. 



Note that the Schwarz inequality implies that D > 0. 
We are in the position now to present the the basic BSA 
theorem: 

Theorem 2 Given the set V of product vectors 
\^a, fa) G R{p), the matrix p* — J^a^aPa is the op- 
timal separable approximation ( OS A ) of p if 

• all Aa are maximal with respect to pa = p — 
J2a'^a^a'Pa', and to the projector Pa; 

• all pairs (Aq,,A^) are maximal with respect to 
Pap ^ P - Y.a'^a,i3^a'Pa' , and to thc projection 
operators {Pa,Pp). 

If V is the set of all product vectors in R{p) (in gen- 
eral uncountable) then the same theorem holds for the 
BSA (for the detailed proof see Appendix to Rcf. po|). 
All information about entanglement is included in the 
matrix 6p. If Sp does not vanish, i.e. if p is not sepa- 
rable, the range R(5p) cannot contain any product vec- 
tor. The reason is that one can use projectors on prod- 
uct vectors that belong to R{5p) in order to increase A. 
The rank of the matrix 5p must be smaller, or equal to 
(M — 1)(A^ — 1). This is because the set of all product 
vectors in the Hilbcrt space H of dimension M x N spans 
a {N -\- M — l)-dimensional manifold, which generically 
has a non- vanishing intersection with linear subspaces of 
H of dimension larger than (TV — 1) x (M— 1). In fact, we 
have proven rigorously that this is the case for 2 x iV sys- 
tems in Ref. ]l5[ | , and presented some rigorous arguments 
for the case M x is Ref. 

In particular, for the case of M = A^ = 2, (5p is a 
simple projector onto an entangled state. For the 2 qubit 
systems it is easy to prove that the BSA decomposition 
is unique and has the form: 



p ^ Aps + {1 - A)P,; Ae[0,l] 



(6) 



where ps is the normalized density matrix. If it had not 
been so, we could have another BSA expansion, lets say 
p — Aps -f- (1 — A)Pe. But taking the convex combi- 
nation of these two decompositions, we obtain another 
BSA decomposition with the remainder dp being given 
by a convex combination of Pe and Pg- Such remainder 
would have then rank 2, and would necessarily contain 
product vectors in its range ||^. If this happened, we 
would be then able to increase the BSA parameter A by 
subtracting from Sp projectors on product vectors in its 
range. That is, however, impossible since A is already 
maximal. For the case of arbitrary dimensions the OSA 
and BSA decompositions are also unique. We present the 
proof of this fact in Section V of this paper. 



3 



III. THE BSA REMINDER OF ® QUANTUM 
SYSTEMS: IS IT MAXIMALLY ENTANGLED? 

We have seen that the BSA reminder of (S) quan- 
tum systems is just given by a projector onto a entan- 
gled state I'f/'e)- This fact is essential and allows to obtain 
the BSA decomposition for some states analytically . 
For many families of states considered by Englert and 
his collaborators the BSA remainder consists of a max- 
imally entangled state. Similar conclusions follow from 
the numerical analysis of Ref. In this section we 

ask therefore a natural question: under which conditions 
the BSA remainder is, or is not maximally entangled? 
Strictly speaking we present here a necessary condition, 
that the BSA decomposition for a generic density matrix 
must fulfill so that the BSA remainder is not maximally 
entangled. 

We concentrate here on generic quantum states which 
have the maximal dimension of the range (r{p) = 
r{p^^) = 4). Let us assume that the density matrix p 
has the BSA decomposition 

p = Aps + {l-A)P^,, (7) 

so that its partial transposition with respect to Alice's 
system, p'^-* = Apf^ + (1 - ^)P^'^ ■ When A is not equal 
to 1, p is entangled, and p'^^ must not be positive definite. 
Let us first observe 

Lemma 3 If p acting in ®C'^has the BSA decompo- 
sition p Kps + (1 - A)P^^, then r{p'^^ ) < 3. 

Proof: Had the range of pj^ been full, one could always 
replace 1 — Aby (1 — A— e), keeping Ap'^^ +eP^^ positive 
definite, while p'^ — ps + eP^_ separable. □ 

The fact that the rank of is not full implies that 
3|i>), such that p^^\v) = 0. Since Pj'* has 3 positive 
and one negative eigenvalue where the eigenvector 
corresponding to a negative eigenvalue in a conveniently 
chosen basis can be written as 



/ \ 

1 

-1 

V J 



, then 7^ 0. If it was not the case, one could also 

replace 1 — Aby(l — A — e), keeping ApJ-^ positive. 

Let us now discuss the optimization procedure, that 
sometimes allow to increase A in the decomposition (|^). 
A given decomposition of such a form is optimal, if it can- 
not be optimized. It will turn out that the optimization 
strategy works only provided "06 is not maximally entan- 
gled. The necessary condition, that the BSA remainder 
is not maximally entangled, is that the decomposition 
cannot be optimized in the sense formulated below. Our 



aim is to formulate this necessary condition in an explicit 
form in this section. 

Optimization procedure: Let us observe that we can 
always write 

|^e)=iVl|ei,/l)+iV2|e2,/2) 

, for any basis |ei),|e2), where (ei|ei) — (e2|e2) — 1, but 
(ei|e2) does not have to be zero. Let |ei), |e2) denote the 
basis biorthogonal to |ei), |e2); we obtain then 

(eilVe) = A^2(ei|e2)|/2) 
(e2|V'e) = A^i(e2|ei)|/i) 

Requiring that — (/2I/2) = 1 the above equations 

allow to determine uniquely A^i,iV2,|/i) and I/2). With- 
out loosing the generality we may assume A^i > N2- Let 
us introduce 



where 

N{af = a^Nl 



-^(aiVi|ei,/i) + liV2|e2,/2)), 



^iV|-t-27ViA^2i?e((ei|e2)(/i|/2)). 



We can now rewrite the BSA projector 

(8) 

We would like to replace the projector P^^ by the expres- 
sion (H) and in this way improve the BSA decomposition. 
To this aim we require that N{af' < 1 which implies that 

a'^Nf + -^N^ < Nl + iV|. Defining now a; = we see 

that N{aY < 1 provided x < < 1. That is only pos- 
sible if Nl ^ N2- The latter conditions fulfilled if tpe is 
not maximally entangled, as described in the following 
lemma: 

Lemma 4 // |?/'e) = iVi|ei,/i) -I- iV2|e2,/2), where 
(ei|ei) = (e2|e2) = 1, then Ni = N2 if ipe is maximally 
entangled. 

Proof: Let us consider a basis in which ji/^e) — a|00) + 
Vl — and assume a general form of |ei) = 

In the basis considered 



we can easy calculate that 



(eilV'e) = a^^lO) + Vl-«2v/T^|l)e-''^, 
so that 

^ll(ei|e2)P = a2p+(l-a2)(l-p) 



aV + (l-a')(l-p') 



(9) 
(10) 



(11) 
(12) 
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e2|ei)P, so that indeed A''^ 



2 11 u, — 2^ that is when the state |V'e) is maximaUy 



Note that |(ei|e2)P 
7V| if a2 = \ 
entangled. □ 

Now we can easily prove 

Lemma 5 If p has the BSA decomposition (Q), then ei- 
ther ijje is maximally entangled, or r{pa) = 3 

Proof: Suppose that r{ps) = 3. If is not maximally 
entangled, the optimization procedure allows to optimize 
the decomposition by taking < 1, but very close to 
one. We can indeed improve BSA for p, provided we can 
subtract -Pe;/2 from Ap^-'^. This means that |e2,/2) 
must belong to the range R{p^^). That in turn requires 
that if \v) = le^, hi) + |e2, we then need (/11I/2) — 0, 
or in another words 



0. 



(13) 



It is easy to see that this equation has many solu- 
tions: for example take |e2) = |ei) and |ei) proportional 
to (^) — |0) + a\l), then the above equation implies 
that"[(w|0) + a* {v\l)][{0\'ipe) + a* (IjV'e)] = 0, which is a 
quadratic equation for a* which obviously has solutions 
for |e2) 7^ |ei). We conclude that either r{ps) = 3, or 
Ni = N2. The latter can occur if and only if is fully 
entangled. □ 

Therefore we have to consider the case r{ps) = 
r{p\^) = 3. From the results presented in the Appendix 
A we know that there exists such a one dimensional fam- 
ily of product states |e2((5), /2(5)), where S is real, such 
that |e2(<5),/2(<5)) e R{ps) and |e^(<5), /2(<5)) € is 
satisfied. 

Now we are in the situation where we can explicitly 
check whether the vector lipe) in the BSA remainder can 
be non maximally entangled. If {ipe) is given and we have 
|e2, 72) = |e(<5), f{S)) for a given ps, then we can calculate 
1/1) and |ei) by 



1/1 



ei 



(e2|V'e) 
l(e2|V'e>|' 

l(/2|V'e)|' 



(14) 
(15) 



and from = 1, we obtain lA^^ij = i^^g^it^^j • Since we 

know now |ei), we can also easily calculate |A'2| = 

\{ei\e^)\ ■ 

We see that the coefficient iVi and N2 can be explicitly 
constructed from ps and \ipe)- We obtain therefore the 
main result of this section 

Theorem 3 If a generic {r{p) = r{p'^^) = 4) state 
rho in (g) has the BSA decomposition p = Aps + 
(1 — A)P^^_, then either ip^ is maximally entangled, or 
r{ps) — r{p'^'^) — 3, and for any expansion of |'(/'e) = 
A^i|ei,/i) + A^i|e2,/2), such that |e2,/2) e RiPs) and 
|e2, /2) e i?(p^'*) holds, it must follow that Ni < N2. 



Proof: The proof is obvious using the lemmas of this 
section, and the optimization procedure. If there exist 
|e2(^), /2(<J)) such that A^i > A'^2, the optimization pro- 
cedure can be applied, which contradicts the optimality 
of the BSA. □ 



IV. ENTANGLEMENT MEASURES 

Before we turn to the main results of this paper let us 
also remind the reader in this section some basic facts 
about entanglement measures and their properties. 

Once one has the physical picture of entanglement as 
a resource, one needs to formulate this concept mathe- 
matically. One way leads through a definition of non- 
entangled, i.e. separable states as discussed in previ- 
ous sections. Another possibility is to try to quantify 
amount of entanglement for a given mixed state. The 
latter approach is realized by defining entanglement mea- 
sures p3t , and by specifying physical properties which 
the entanglement measure should have. There are sev- 
eral versions of definitions of the entanglement measures; 
here we follow the approach of Plenio and Verdal : 

Definition 3 Let p be a quantum state acting in a 
Hilbert space Hab = 'Ha®'Hb, then the function E{p) 1— > 
R is called entanglement measure if it satisfies: 

1. E{p) — Q, if p is separable; 

2. Local unitary operation leave E{p) invariant, i.e. 
E{p) = E{Ua ® UbpU\ ® Ul); 

3. Let AiA^^ ® BiSj = 1 be some complete lo- 
cal measurement (i.e. local positive operator valued 
map (POVM)), then 



E{p)>Y,T^r{pi)E{pjTT{p,)), 



(16) 



where pi :— Ai(^ BipA\ ® Bj . This property means 
that entanglement measure cannot increase in the 
mean under local operations. 

4-. For pure states the measure of entanglement should 
reduce to the entropy of entanglement, which 
is defined as von Neuman entropy of the reduced 
density matrix, pA = Tt^bP (or, alternatively pb = 
TrApj, 



E{p) := -Tr(pAlnpA); 



(17) 



5. Entanglement measure should be additive which 
means that 



E{pi^P2)^E{pi) + E{p2). 



(18) 
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It should be pointed out that the necessity of the last 
two conditions is still disputed in the literature p5| , p6| , 
and therefore we will just concentrate on the first three 
conditions. Notice, that in Eq. ( p^ it may happens that 
E{p,/lY{p,)) < E{p). 

To complete this section, let us list some of the most 
widely used entanglement measures. Typically, they ful- 
fill some, but not all of the conditions 1-5 of the Def. 



1. Entanglement of formation is defined as 
EF:^uimY^p,S{p\), (19) 



where S{pa) — Tr(pyilnp^) is the von Neumann 
entropy and the minimum is taken over all the 
possible realizations of the state, p = J2i 
where p\ — TrB(|?/)j)(?/)j|). Notice that in the case 
where p is a pure state [p = , the von Neu- 

mann entropy of the reduced density matrix is an 
entanglement measure. The physical meaning of 
the formation measure is the minimal amount of 
pure state entanglement needed to create a the 
given entangled state. Calculation of Ep for a given 
state is a very difficult task. Remarkably, Wooters, 
has derived the analytic formula for Ep for an ar- 
bitrary two qubit state [p7| . 

2. Relative entropy entanglement measure ||2j] 

is defined as 



min E{p\\ps); 

Ps 



(20) 



where the minimum is taken over all separable 
states Ps and E{p\\ps) is the relative entropy, which 
is given by the expression 



E{p\\ps) :=Tr(p(lnp-lnp,)) 



(21) 



3. Bures entanglement measure is defined as 
E{p) :=min(2-2v/F(p,p,)), (22) 

Ps 

where F(p, ps) is the Uhlmann's fidelity F{p, ps) := 
(Tr(yY^pp7\/p))^- This entanglement measure 
does not fulfill the last two conditions of Definition 
3. 

In the recent years a very promising approach has been 
initiated by Vidal who has shown that more parameters 
(the so called entanglement monotones) are required in 
order to quantify completely the non-local character of 
bipartite pure states [p6|. 



V. THE BSA ENTANGLEMENT 

Let us now investigate how do the local POVM's influ- 
ence a given BSA decomposition. To this aim we consider 
a POVM of the form of '£^ V^Vt =1. ^ Ai ^ Bi. 
After the i-th result is obtained in the measurement we 
obtain the following density matrix 



Pi 



Tt{V,pV^) 



A 



Tr(X^,pl//)'? Tr(T/,p,yJ) Tr(l/,P„y/ 
(1-A?^)(^)- 



Defining now 



A, := A 



5pi ■-- 



Tt{V,PsK^) 

Mv.pv,^) ' 

Tr(y,P^V-;) 
" Tr(y,p,^/) ' 



Tr(K:Pal^/ 
V,6pV^^ 



MVrSpV.^y 
We rewrite the result as: 

V^pV,^ ^ Pi^ A, J2 ^iaP^c + (1 - A,)(5p, 

Q 

We observe that 

l-A = ^(l-A,Tr(y.py;)) 



(23) 



holds. Since for the BSA decomposition of pi the inequal- 
ity 

AssA. > A, (24) 
holds, we get from (^Sj) that 

l-A>^(l-AssA.Tr(l/,pl/;)). (25) 

i 

The result (^) allows to prove the following property: 
Property 1 The BSA entanglement measure 

E{p) = l~ABSAip) 

fulfills the properties 1.-3. of the Def. |^. 
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Proof: 

1. If p is separable, i.e. p ^ ps then A = 1, and 
E{p) = 1 - A = 0. 

2. U p = Ua® UbpU\ ® then obviously E{p) > 
1— A = E{p), and vice versa, since we can invert Ua®Ub- 
That means that E{p) is invariant with respect to local 
unitary transformations. 

3. Finally, if we apply a local POVM, we obtain 

E{p) = 1 - A > ^(1 - ABs^,Tr(F,pF;)) 

i 

>Y,E{p,)Tr{V.pV^), 



where p^ = V,pV^ /Tt{V,pV^). This follows from 

It is worth noticing that the above argument holds for 
the Hilbert spaces Ha ® Ti-B of arbitrary dimensions. 



VI. THE UNIQUENESS OF THE BSA 

In this Section we turn back to the general case and 
present a proof that the BSA in any Hilbert space is 
unique. To this aim we prove first a lemma, and that the 
major result. 

Lemma 6 Let a hermitian density matrix p has a de- 
composition of the from p — Aps + (1 — A)Sp, where 
Ps is the separable part which has the structure ps — 
A^^^j^ AqPci, with Pa being the projection operators 
onto the product states \ea, fa) and X]q=i = 1- Then 
the set of {A^}, which are maximal with respect to the 
density matrix p and the set of the projection operators 
{Pa}, form a manifold which generically has a dimension 
n — 1 and is determined by the following equation 

n n n 

1 - J2 ^^Di + XI ^i^'jDi] - A,Aj AfcDyfe + . . . 

i i<j i<j<k 

+ Ai,A,, ...Ai„,Al^2...^™ + 

il<i2<-<lm 

... +(-)"AiA2...A„i?i2...„-0 (26) 

where the set o/ {-Di^ia - i™ } '^'^^ subdeterminants (mi- 
nors) of the matrix D, which is defined as 



D = 



{i>2\p-^\i^l) (^2|p-l|^2> 



(V'llp-^IV'n) \ 

mp-'Hn) 



\mp-^i^l) {Mp-^4>2) ... (^„|p-^|V'n) / 



and where by {\'ipi)} we denote for shortness the product 
vectors which are building the projection operators {Pi = 



Proof: Let us first remark that generically the matrix D 
does not have a block structure. If the matrix D consists 
k diagonal rife-dimensional blocks, then not only Eq. ( [2^ ) 
is fulfilled, but also the k corresponding equations for 
the blocks, so that the corresponding manifold has the 
dimension rife , and is a cartesian product of k manifolds 
of dimension — 1. In the following we will concentrate 
on the generic case. 

The proof of the lemma goes with induction. First we 
prove it for n = 2 and we get 

1 - AiDi - A2D2 + AiAzI^ia = 0, 

or for n = 3 where we get 



1 -AiDi - A2D2 - A3D3 + A1A2Z312 
+ A2A3D23 - A1A2A3D123 = 0. 



AiA3i? 
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Now, let us assume that the lemma is true for n, and 
show that it must also be true for n + 1 . Let p has the 
decomposition p = Aps + (1 — A)Sp, with 



ri+l 



Ps = A^ AaPa. 



The lemma holds for the matrix p — p ~ 
A„+i|'0„+i)(V'n+i| so that the first n coefficient A^ fulfill 
Eq. ( p6| ) with coefficients D calculated as above with the 
substitution p^^ ^ p^^ = (p - A„+i|V'„+i)(^/'n+i|)"^. 
The latter inverse can be calculated using power series 
expansion in the projector A„+i|'0„_).i)(?/',i+i|. The re- 
sult is 



(P-A 



i|V„+i)(V^„+i|)~V.) = P^V.) + 



A„+l \p~^\i)^){i)^\p~^\'\l}n 



1 - A„+i (V'n+l \p ^ IV'n+l) 

Inserting the above result to equations defining the sur- 
face for the first n A's we get, after tedious, but elemen- 
tary algebraic calculation 



1 -^A.A 



^A^AjA,, - A.AjAkD. 



E 

n <i2 <---<lm 



A^^Ai^ 



+ (-) 

n< 

+ (-)" E 

il<i2<-- 

+ (-)"+iAiA2...A„+ii? 



i<j<k 



.A,;„A 



ijk 



+ 1-L^12...n+ 



1=0 



which proofs the lemma for n -|- 1 . □ 

Note that in particular, if the decomposition discussed 
in the above lemma is the BSA, then the correspond- 
ing A's fulfill Eq. (^6|). This observation allows us to 
prove the uniqueness of the BSA in arbitrary dimension. 
It is important to note that the surface defined by Eq. 
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( Pq ) can be considered for arbitrary A's, not necessarily 
positive! This surface is strictly convex and divides the 
space of all A's into two sets: a convex set of those sets 
of {A's} which have the property that p — A^^jt^ AqPq 
is positive definite, and concave set for which the latter 
matrix is not positive definite. If this surface contains a 
part of a hyperplane (linear subspace), it must contain 
the whole hyperplane, since it is defined by the polyno- 
mial equation (^. This observation is essential to prove 
the uniqueness of the expansion. 

Lemma 7 (The uniqueness of the BSA) Any den- 
sity matrix p has a unique decomposition p — Ap^ + (1 — 
A)Sp, where ps is a separable density matrix, Sp is a in- 
separable matrix with no product vectors in its range, and 
A is maximal. 

Proof: The proof the lemma goes by assuming the 
decomposition is not unique; then there must exist at 
least two BSA decompositions, p = Apsi + (1 — A)Spi 
and p = Aps2 + (1 ~ A)Sp2, with the same maximal A. 
Now, any convex combination of these two BSA decom- 
positions is also the BSA decomposition, 

p = epsi - (1 - e)ps2 + eSpi + (1 - e)Sp2 
= ^(eAAn - (1 - e)AA2,)P» + [eSpi - (1 - e)Sp2) 

i 

= Ps{.f) + 5p{^), 

where e S [0, 1]. The part of the one dimensional hyper 
plane (line) eAii — (1 — e)h.2i for e G [0, 1] lies on the 
surface (26). 

From the form the surface it follows that the whole line 
eAii — (1 — e)A2i for all e lies on that surface. This cannot 
be, since for some e ^ [0, 1], and 5pi ^ 6p2, <5p(e) must 
become nonpositive definite. This is easy to see since for 
e ±oo, (Sp(e) (X bp\ — 5p2, and the latter matrix is 
non zero and has the trace zero, so that it has to have 
eigenvalues of opposite signs. This is thus a contradiction 
with the assumption made at the beginning, ergo the 
BSA decomposition must be unique. 



VII. THE PPT BSA 

In this section we discuss in detail generalization of the 
BSA approach for PPT states used in Refs 

Theorem 4 Let p he a arbitrary PPT state. For any 
countable set V = {Pi = l^i, /i)(ei, such that 
|eij/i) S R{p) '^''^d. \e*,fi) E R{p^'^), there exists the 
best separable approximation of p in the form 



A 



PPT 



maxv(Tr(p,[F])). 



(28) 



p^Aps + {l- A)Sp, 



(27) 



where ps — AiPi is a separable state, A is maximal, 
and both dp > 0, and Sp^-"^ > 0. We call such a de- 
composition a PPT BSA if it preserves the PPT of the 
remainder Sp and 



Proof: Let us consider the set of all separable matrices 
Ps = Si'^i|ei,/i)(ei,/i|, where \ei, fi) G V,p-ps > and 
p-^'* — pj''^ > 0. This set of p's form a convex and bounded 
set, which means that this set is compact. Because of 
the compactness there must exist a separable matrix ps 
which has maximal trace A = Tt{ps[V]). By expanding 
V we will finally get the maximal PPT contribution. □ 

Let us analyze the PPT BSA decomposition in more 
detail. All information about the PPT entanglement is 
included in the PPT BSA parameter A and 5p. If the 
PPT BSA remainder 5p does not vanish, then there exists 
no product vector |e,/), such that |e, /) S R{5p) and 
simultaneously |e*, /) S R{6p^^) is satisfied. This means 
that the PPT state dp is entangled. 

We introduce now, just like in the first version of the 
BSA, a procedure of constructing the matrix p^. But 
before we do this let us define some basic concepts for 
that: 

Definition 4 A non-negative parameter A is called 
PPT maximal with respect to a positive PPT oper- 
ator p, and a projection operator P — G ^ */ 
p - AP > 0,p'^'^ - Ap^-* > 0, and for every e > 0, the 
matrix p — {A + e)P is not a PPT state. 

This means that the A is the maximal contribution of 
P that can be subtracted from p by maintaining the 
PPT of the difference. Now let us introduce the following 
Lemma: 

Lemma 8 A is PPT maximal with respect to p and P ~ 
|e,/)(e,/| iff: 

. if |e,/) ^ R{p) and |e*,/) ^ R{p^-), or |e,/) ^ 
R{p) and \e*, f) £ i?(p^^) or |e,/) S R{p) and 
|e*,/) ^R{p^'^) then A = 0; 

• if\e,f) e R{p) and |e*,/) G R{p^^') then 



A = min(((e,/|i|e, /))-!, ((e*,/|-i^|e*,/))-i 



(29) 



Proof: 



From lemma 

-1 



(|^) we know that A 



((e, /|i|e, /))"^ is the maximal contribution to p and 
A — {{e* , f\jrx\e* , f))~^ is the maximal contribution to 
P"^-*. In order to maximize and keep the PPT of the 
difference we have to take the minimum of A and A.n 

Definition 5 A pair of non-negative (Ai,A2) is called 
maximal with respect to p and a pair of projection op- 
erators Pi = |ei,/i)(ei,/i| and P2 = |e2, /2)(e2, /2I iff 

• p- AiPi - A2P2 > and (p - AiPi - A2P2)*'* > 0, 
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• Ai is PPT maximal with respect to p — A2P2, 

• A2 is PPT maximal with respect to p — AiPi, and 

• Ai + A2 is maximal. 

The conditions for PPT maximizing of pairs Pi — 
|ei,/i)(ei,/i| and P2 = |e2, /2)(e2, /2I are described in 
appendix ^. 

Let us now prove that for a given countable set V of 
product vectors we can obtain the optimal PPT separable 
approximation by maximizing all pairs of productvectors 
in V. But before we do this, we have to define the PPT 
BSA manifold: 

Definition 6 Let the equation F{Xi, . . . , \k) = (or 
Ai — /i(A2, . . . , Afe)J describes the BSA manifold with 
respect to p, and F(Ai, . . . , A^) — (or Ai = 
/i (A2 , . . . ,Xk)) for p*-* . Without loosing generality in or- 
der to obtain the manifold which preserves the PPT of 
the differenz {p — ps) we have to define 



^(A. + /,)U=A 




x^A < 0, (32) 



Ai = min (^Ai = /i(A2, . . . , \k),Xi = /i(A2, ■ ■ ■ , >^k 
= /i(A2,...,Ak). (30) 

The implicit form will then be given by F{Xi, . . . , Xk) = 
0. 

Notice that the PPT BSA manifold is contineous and all 
most everywhere differentiable. 

Tiieorem 5 Given the set V of product vectors \ei, fi) G 
R{p) where also \e* , fi) G R{p'^'^), then the matrix ps = 
^j^-^ AiPi is the optimal PPT separable approximation 
of P «/•■ 

• all Ai are PPT maximal with respect to pi = p — 
'^i/^i Ail Pi' , and to the projector Pi; 

• all pairs (Ai,Aj) are PPT maximal with respect 
to Pij = p — j -^i' -Pi' ' '^'^'^ to the projectors 

Proof:If Ps is a PPT BSA decomposition then all A^, 
as well as all pairs (A^, Aj) must be PPT maximal (oth- 
erwise maximize Aj would increase the trace of ps). 

To prove the inverse, consider matrices ps = ^iPi 
for which all individual Ai are PPT maximal. This means 
that Ps belongs to the boundary of the set Z of all sepa- 
rable matrices such that p — Ps > and {p — ps)*^ > 0. 
This boundary is the PPT BSA manifold: 



F(Ai,...,A/f) = 0. 



(31) 



The manifold ( pT| ) can be written as a function A^ = 
fii{^j}j¥=i)^ depending on which size of the manifold we 
are. Let = J^i AiPi be the separable matrix for which 
all pairs of A's are PPT maximal. The maximum of 
(Ai, Aj) then implies that 



for all sides of the manifold F = and i,j. This means 
that p™ is either a local maximum or a saddle point (not 
necessary the same derivative in every direction of A = 
A). Now we have the same situation just like in the 
original version of the BSA. The later possibility cannot 
occur, since the set Z is convex (i.e. if ps,p^g Z then 
eps + {l~e)p's e Z for every < e < 1). Since p2| describes 
also a convex set it can for sure not be a saddle point. 
The same argument holds also for the local minimum. 
And finally the local maximum must be also a global 
one, because on a convex set there can not exists two of 
them. This means that ps = p™.a 

One should stress out at the end of this section, that 
the PPT BSA can be straight forward generalize to mul- 
ticomposite systems. 



VIII. CONCLUSIONS 

In this paper we have presented several novel results 
concerning the BSA decmpositions of density matrices 
of composite quantum systems. General results concern 
the uniqueness of the BSA decompositions, the existence 
of the BSA entnaglement mass, and the efficient meth- 
ods of construction of the BSA decomposition for PPT 
states. More specific results for two qubit systems deal 
with the necessary conditions, that the projector onto 
a nonmaximally entnagled state proviedes the remain- 
der in the BSA decomposition. There are several open 
questions concerning the BSA decompositions in higher 
dimensional Hilbert spaces: what is the structure of re- 
mainder in such a case, how to parametrize the remain- 
ders (tha so called edge states 18 1 in the case of PPT 
BSA). The physical interpretation of the BSA entangle- 
ment mass is not known so far. In the case of 2 x 2 space, 
we hope that our results, togehter with remarkable ana- 
lytic results of Englert and his colleagues |^l| will bring 
us closer to the challenging goal of analytic construction 
of the BSA decomposition for arbitrary two quibit den- 
sity matrix. 

This work has been supported by 
Deutsche Forschungsgemeinschaft (SFB 407 and Schwer- 
punkt "Quanteninformationsverarbeitung"), and by the 
1ST Programm "EQUIP" . We thank D. Bruss, J.I. Cirac, 
B.-G. Englert, P. Horodecki, B. Kraus, A. Sanpera, R. 
Werner and M. Wilkens for fruitful discussions. 



APPENDIX A: PRODUCT VECTORS IN THE 
RANGE 

In this appendix we prove some lemmas that has 
been used in the section IV. Both the results as well as 
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the proofs are very much parallel to the one used by 
Woronowicz [lO| . 

Lemma 9 If p is a density matrix in a 2 x 2 space hav- 
ing a positive partial transpose and r{p) = r{p^^) — 3, 
then there exist a product vector |e, /) G R{p) such that 
\e*J) eR{p^-). 

Proof: Let there be given a density matrix p — 
A B \ 

^1 ^ j (A and C are invertible, because otherwise 

we would have a product vectors in the kernel [ p5[ , 
and the existence of |e, /) would follow from the re- 
sults of Ref. ||l5|). Now, we choose the basis in Ha to 

tain that B{a*) 



)}. In this new basis we ob- 



:(1 



A B 
C 



a function of a* only. This means that we can choose a 
such that deti3(a*) — deti3^(a) = 0. Choosing such an 
a, we get r{B) = t{B*) — 1. 



The next step is to perform a non unitary, but invert- 
ible local transformation p ^ Ia® '^P-^^ ® ^^'^ 
redefine A ^A^, B 



After that, the 



new matrix is given by p = 

our assumption that r(p) — 3, 
A 



Now, we use 



1 n 1 

Vc Vc- 
A B 
St / 

from which it follows that 
BB"^ + AP, where P is a projector on some vector 
\ip). The assumption that also r{p^^) — 3, leads us to 
A = B^B+XP, where P is a projector on some other vec- 
tor IVi). This leads us to BB^+XP ^B^B+XP, and since 
tri^BB^ -B^ B) = 0, we get that A = A. What is the neces- 
sary condition now for (J^^^) G ^{p) and (^Ij]^^) € r{p^^) 
? This condition means nothing else than that there exist 
two vectors, lets say 



and (j^j) , such that 



PPt + XP B 
Pt / 

B^B + XP Pt 
P / 



l/> 
I/) 



from which we get the equation 



1 



1 - zP 



\lp) = T] 



1 



1 -z*Pt 



(Al) 
(A2) 

(A3) 



with some complex rj. In order to proof our lemma we 
must show that there exist a solution for ([A^). The trick 
is now to describe the right side of the equation (A3) 
as a complex conjugate of the left side, so that we can 
construct a solution explicitly. 



We will show now that the equation (A3) can indeed be 
transformed into the form 



1 



1 - zP 



1 



1 - z*B 



(A4) 



where is the Pauli matrix. Defining = ("^), 

This 



we must have that vi = rie'^'''v2 and V2 
equation has a solution if vi — ve^^ and V2 — ve^^^^ , 
where = ||w2|| = v. Lets take now an arbitrary 6 
and require (^i^) ^ I''/')' which means that 



js 



1 



1~ zB' 







(A5) 



must hold. Obviously, this equation has not only one 
solution, but an infinite family of solutions for every S. 



Let us now proof that equation (A4) indeed holds. 
First we choose a basis \ipi) ,\'4'2) such that B^ B — 

^ ^ Therefore we have that X{P ~ P) = 



BB^ 

A 




-A 





-A 



Since the overall phases of 1-0) and \tp) are 



irrelevant, we parameterize lip) and |0) in our new basis 
as IV') = (yr^e-*): IV') = (^^)- This parameteriza- 
tion yields p = p, (p = (p and A = A(l — 2p). We ob- 
serve now that there exist always a unitary K such that 
KBK^ = B'^. From this trivially follows of course that 
(K'^YB^K^ = P, and therefore {K^)^KBK^K'^ = P, 
from which then BU — UB, where U = K^K^ . 

Now we will proof that K = e**^" {^10 

Let M = BB^ - B^B = A(P - P) (Note that M = M* 
in our basis). Then we have KMK^ = P^P* -B*B^ = 
B*{B'^)* - {B^)*B* = -M* = -M. Therefore M = 
XiKliPXiPlK^ - K\'i/j){i:\K^), and for the vectors |V'),|V') 
we get 



Km 



e'f'y/pe"'' 



This implies K — 

(j),9i + tp = ipi,ip2 = 

if 01 02 then U = 







and therefore O2 ^ fi + 
= 02. But, 
U will 



and 939 + 



e-H«i-«2 
commute with P, if P is diagonal in the chosen basis. But 
then PPt - Ptp = 0, from which follows that \^) {ip), 
and thus (p) in the range of p which proves the Lemma. 
This means that di —62, and K ~ e^^^ax- Since the 
overall phases of K are irrelevant, we can assume that 
K — Gx- This proves however (A4), which consequently 
proves the Lemma too. 

The reader made think now that we have finished the 
proof of the Lemma, but remember that at the beginning 
of the proof we have made a non unitary local operation. 
What we must do now is to retransform the density ma- 
trix p, and check if our results after that still holds. Let 
us see what happens after the inverse transformation: 
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_ / VcBB^Vc + xVcpVc VcbVc \ 

VCB^VC C ) 

Demanding that (jj'^^^) e R[p) and (^'j^^) e R{p^^) leads 
to the foUowing conditions: 

T^TT^^I'/-') = ^1 1 

- > = \/c7?(i - r{z)B^m, 

{I- f{z)Bm=r^{l- r{z)B^)<TA^)- 

We see that the equations are equivalent after the rescal- 
ing, so that the Lemma holds . □ 

The prove of the above Lemma allows to parameterize 
the set of all product vectors \e{5),f{5)), which satis- 
fied the condition \e{5),f{5)) G R{ps) and \e{5)*,f{5)) £ 
R{p'^'^), by an one dimensional real parameter 6. This 
will be used in Section IIL 

APPENDIX B: PPT PAIR MAXIMIZING 

In this appendix we explain how to PPT max- 
imize a pair of product projectors {\ipi){ipi\ = 
|ei,/i)(ei,/i|,|MV'2| = |ei,/i)(ei,/i|). 

As we know from the BSA, the BSA manifold for p 
and (l^iXV-il = |ei,/i)(ei,/i|,|^2)(V'2| = |ei, /i)(ei, /i|) 
is given by 

F(Ai, A2) = 1 - AiDl - A2D^ - AiA2D° = 0, (Bl) 

where - (ei, /i|p-i|ei, /i),Z?§ - (ea, Mp-^lea, /a) 
and = (ei,/i|p-i|ei,/i)(e2,/2|p-i|e2,/2) - 

^|e2,/2)|P- But also we have to consider the 
BSA manifold for p^^ . This one is given by 

F(Ai , A2) = 1 - AiDl - A2DI - Ai AaZJi = 0, (B2) 

where 

Dl = {elJi\ip'-r^\el,h),Dl = (e*, /2|(p*-)-^|e^, /2> 
and = {elJi\i/-)-'\el,h){e*2j2\{p'-)-'\e*2j2)- 
^|e27/2)|P- Now we have to consider two 
basic cases which can occur. 

Case 1: One of the BSA manifolds is under the other 
manifold. Without loosing generality we assume that this 
is = 0. Then we have the situation just like in figure 
1. 




In that case we have to take the maximum on the man- 
ifold F = 0. From lemma ^ we know the condition for 
that. Of course we are also including in the case 1 that 
there can be an overlap at one endpoints (i.e. ii = 

Case 2: The manifolds have a cross section point be- 
tween < Ai < max ^tjd-j tjt^ • Without loosing gener- 
ality we assume that this describes Figure 2. Now we can 
see from Figure 2 how the PPT BSA manifold = is 
constructed, and why it is not differentiable every where. 




FIG. 2. The manifolds have a cross section point As 

Let us denote by A„j the maxima of the manifold F = 
and also Am as the maxima of = 0. Now we can have 
the following situations: 

• If Am < As and Am < As then one has to take 

Aniax — Am-, 

• If A,„ > As and A„i > As then one has to take 

Amax — Am-, 

• If Am > As and Am < As then one has to take 

Amax — 5 

• Both maxima are in As, so that Amax = As- 

• The case where Am < As and A„j > As can not 
occur; 
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